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For the quantum depolarizing channel with any finite dimension, we compare three schemes for 
channel identification: unentangled probes, probes maximally entangled with an external ancilla, 
and maximally entangled probe pairs. This comparison includes cases where the ancilla is itself depo- 
larizing and where the probe is circulated back through the channel before measurement. Compared 
on the basis of (quantum Fisher) information gained per channel use, we find broadly that entan- 
glement with an ancilla dominates the other two schemes, but only if entanglement is cheap relative 
to the cost per channel use and only if the external ancilla is well shielded from depolarization. We 
arrive at these results by a relatively simple analytical means. A separate, more complicated analysis 
for partially entangled probes shows for the qudit depolarizing channel that any amount of probe 
entanglement is advantageous and that the greatest advantage comes with maximal entanglement. 



PACS numbers: 03.67.-a 



I. INTRODUCTION 



A quantum operation, or channel, T is a trace- 
preserving, completely positive map for representing the 
change of state of a quantum system. Quantum channels 
are the fundamental building blocks of quantum informa- 
tion processing. In any physical implementation, though, 
the channel is usually not fully known and must be de- 
termined experimentally by observing its effect on pre- 
pared quantum systems, or probes. This is the problem 
of quantum channel identification [H-01j a l so known as 
quantum process tomography , and a major challenge 
is to determine the optimal quantum state and measure- 
ment for the channel probe. In particular, the probe 
may be entangled with an ancilla system, or with other 
probes. 

Quantum channel identification is statistical by nature 
and is generally formulated as a parameter estimation 
problem: the unknown channel Tg is supposed to be- 
long to a parametric family {Tg, 9 £ 0} of channels, 
and we identify the channel by estimating 9. The gen- 
eral scheme for doing this is to prepare the channel input 
probe in a chosen quantum state a, make a quantum 
measurement M. of the channel output pg = Tg[a], and 
base the estimation of 9 on the measurement's registered 
result X. This process (see Fig. [IJ is repeated to ob- 
tain n independent, identically distributed measurement 
outcomes Xi, . . . , X n , and then 9 is estimated by an es- 
timator 9 = 9{X\, . . . , X n ). This arrangement, in which 
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a prepared quantum system is used to probe a quantum 
channel, is standard for channel identification, it is com- 
mon in interferometry Q, and it occurs generally when- 
ever a quantum measurement is made indirectly by ob- 
serving a change in an ancillary part of the measurement 
apparatus. 
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6 — 9{X\, . . . , X n ) 



FIG. 1. Set-up for probing a depolarizing channel Fg 
probe system prepared in state a. 



by a 



We study the problem of identifying the qudit (di- 
mension d < oo) depolarizing channel. This channel 
is a standard model of quantum noise, it is analytically 
tractable, and it is the basis of investigation in a range 
of contexts @, Hcil - fl4j . Our study of the qudit depolar- 
izing channel proceeds most directly from the works of 
Fujiwara and Sasaki et al. Fujiwara [l[ considers the 
qubit (d = 2) depolarizing channel with probes prepared 
by three different schemes. Using a likelihood approach 
based around quantum Fisher information, Fujiwara ob- 
tains comprehensive "local/asymptotic" optimality re- 
sults for his three identification schemes, even including 
probes prepared in partially entangled states. Sasaki et 
al. study the more general qudit (d < oo) depolariz- 
ing channel from a Baycsian perspective, comparing the 
same channel identification schemes as Fujiwara. Their 
Bayesian approach yields "global" optimality results for 
any d < oo, though their results for partial entanglement 
are limited to d — 2. We take the same approach as Fu- 
jiwara but treat the general qudit depolarizing channel, 
introducing a common simple expression for the quan- 
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turn Fisher information to address all three identification 
schemes treated in [H, 0] • We extend this method to en- 
tanglement with an ancilla that is itself depolarizing and 
to a scheme we call probe re-circulation, whereby the 
probe is passed multiple times through the channel be- 
fore measurement. By a separate analysis we show for the 
depolarizing channel that any degree of entanglement is 
beneficial for qudits of any dimension and that maximal 
entanglement is most beneficial. This broad usefulness 
of entanglement in any dimension is known in only a few 
instances [l5j . 

The remainder of the paper is arranged as follows. We 
set the stage for our work in Sec. Hi] by recalling the defi- 
nition of quantum Fisher information and its role in the 
quantum Cramer-Rao inequality. We then derive a gen- 
eral expression for the quantum Fisher information in the 
output of the qudit depolarizing channel. In Sec. IIIII we 
apply this expression to channel identification schemes in 
which the qudit probe is unentangled, maximally entan- 
gled with an ancilla qudit, or maximally entangled with 
another probe (the three schemes studied in dHll)- F° r 
each of these schemes, we treat the possibility that the 
probe is circulated back through the channel before mea- 
surement. We further apply our expression for quantum 
Fisher information to a probe that is maximally entan- 
gled with an ancilla qudit that is itself depolarizing. Par- 
tially entangled probes must be addressed with a separate 
method; we do this in Sec. IIVI We take up the question 
of optimal probe measurement in Sec. [V] and summarize 
our results in Sec. IVI1 



II. PRELIMINARIES 

The quantum Fisher information J{9) = J[p$] in the 
parametric output pg = Tg(a) of a quantum channel 
bounds the ultimate precision of the estimation of the 
parameter 9 attainable by quantum measurement Ai of 
pg. The quantum information (Cramer- Rao) inequality 
states for n independent, identically distributed registra- 
tions of any quantum measurement M. and any unbiased 
estimator 9 that 



mse [9] ^ 



1 



iJ{9) 



where mse [9] = E [(6 — 0) 2 ] is the mean squared error of 
9. The quantum Fisher information J{9) in Eq. ([IJ is 



J(9) = J[pg] = tl[p e Ll] 



(2) 



where Lg is the quantum score operator (symmetric log- 
arithmic derivative) of pg defined by 



iLg 



On 



(3) 



classical Fisher information is brought out in [17[ ; and re- 
cent presentations are [TH, [l9j]. The Cramer- Rao bound, 
Eq. JTJ), is asymptotically achievable so the larger J (9) 
is in Eq. ([1]), the more precisely 9 can potentially be es- 
timated. We therefore interpret J(9) to be the overall 
information about 9 available in a parametrically defined 
quantum state pg. In this sense J (9) serves as a quan- 
titative measure of the relative merit of different chan- 
nel identification schemes. The information inequality, 
Eq. Q i an d this interpretation of quantum Fisher infor- 
mation, readily extends [l8[ to biased estimators, depen- 
dent and non-identically distributed registrations, vec- 
tors of channel parameters [20j . and transformations of 
parameter. 

The qudit (d-dimensional) depolarizing channel de- 
fined for any qudit input state a is 



1 - 



9a 



(4) 



where 1 — 9 is the probability of depolarization and 
is the completely mixed qudit state. The channel Tg is 
completely positive for ^ 9 ^ 1 . In fact [2l[ , Tg is com- 
pletely positive for -l/(d 2 - 1) s$ 9 ^ 1; we restrict our 
attention to ^ 9 ^ 1 where 9 reflects solely depolar- 
ization. We take advantage of a general transformation 
of parameter to simply and flexibly express the quantum 
Fisher information J(9) for the qudit depolarizing chan- 
nel. This formulation of J(9) is used in the following 
section to develop and compare schemes for estimating 
9 and identifying Tg. In particular, we generalize results 
known [l[ for the qubit depolarizing channel to the set- 
ting of qudits. 

We will need to consider a variety of dcpolarizing-type 
channels, whose general form is 



v -> ; -fd 



h(9)a 



(5) 



where d is the dimension of the system and $C h(9) ^ 1 
is monotonic and has the interpretation that 1 — h(9) 
is a generalized depolarizing probability. Given that the 
probe state a at the input of the qudit depolarizing chan- 
nel has spectral decomposition a = Y^i=i s i\ s i)( s i\ it fol- 



(1) lows that pg = Tg(a) = £\ =1 Ai(0)K)(si| where 



M0) 



1 - h{9) 



Sih{6). 



(6) 



Only the eigenvalues Aj(0) of pg depend on (h_ its eigen- 
vectors | do not, so pg is quasi- classical pj, HJ. In 
this situation pg and the score operator Lg commute, Lg 
is readily calculated, and we have 



J{9) = Tr [pgL 2 ] = Tr \pg\dgpgf] = 



a,(0) 



(7) 



where dg signifies differentiation. Quantum Fisher infor- 
mation originated with Hclstrom [16}; its relationship to 



Fujiwara [l| observed that the quantum Fisher infor- 
mation J (9) = J[Tg(a)] in a channel output r A1 (cr) about 
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a channel parameter 9 is maximized by a pure state a. 
This observation applies in particular to the qudit depo- 
larizing channel, so we restrict our attention here to pure 
probe states. Let {\si), i = 1, . . . , d} be an orthonormal 
basis relative to which a = \si)(si\. The channel output 
is then 



Pe 



1 + (d - l)h(9) 



si)(si\ + 



1 - h(9) 



5>><*l- (8) 



Reading the eigenvalues of pg from Eq. ([8]), we find that 
the quantum Fisher information Eq. ((7|) in pg is 



J(6) = 



(h'(9)f 



(l-h(0))[h(6) 



i 



(9) 



This expression for J(9) applies to all the qudit depolar- 
izing channel identification schemes of the next section. 



III. PROBING WITH ENTANGLEMENT 

We first consider channel identification in which the 
qudit probe is prepared in an unentangled pure state 
a — we call this scheme O. The probe is applied to the 
depolarizing channel with depolarizing probability 1 — 6; 
i.e., h(9) = 9. We have immediately from Eq. {[9j) that 
the quantum Fisher information available in the channel 
output about 9 is 



Jo(6) 



(1-9) (t 



i 



(10) 



with no input entanglement. This quantum Fisher infor- 
mation is the same for any pure input state <r, because 
the qudit depolarizing channel is unitarily invariant and 
J[pg] = >/[UpeU T ] for any unitary transformation U. 

We now consider an externally entangled probe. 
Specifically, in scheme E one of a maximally entangled 
pair of qudits is passed through Tg with depolarizing 
probability 1 — 9; again h(9) — 9. The state of the max- 
imally entangled qudit pair is a = \p)(p\ with 



Im> = 



sit 



\sdtd) 



Vd 



(11) 



where {|si)} and {\ti}} are orthonormal bases for each 
qubit. The channel output pg — (I ® Tg)(o~) is, after 
some calculation, 



1 - 
~d? 



■L« + 9 a 



(12) 



This output is identically that of a oP-dimensional qudit 
depolarizing channel with depolarizing probability 1 — 9 
and a probe prepared in a pure state. From Eq. (|10[). 
then, the quantum Fisher information in Eq. (|12p about 
9 obtained by a maximally entangled probe is 



MS) = 



i 



(i- 



d^T 



(13) 



and the information gained using a maximally entangled 
probe is 



M6)-Jo{e) 



d{d- i)(i 



(i 



d6)(l 



<poy 



(14) 



This gain is positive for all d < oo and all 9 < 1. 

One might entangle two qudits and use both as probes 
as a scheme for channel identification. In scheme B two 
qudits in the maximally entangled state a = \p)(p\, with 
as in Eq. (fTTj) . are each passed through the qudit 
depolarizing channel. This yields the channel output 
pg = (Tg (g) Tg)(a), in which case 



1 



Pe 



d? 



>-d 2 



9 2 a. 



(15) 



This is identically the parametric family of output states 
that would result from passing a c? 2 -dimensional unen- 
tangled qudit in the pure state a through a channel with 
depolarizing probability 1 — h(9) where h(9) = 9 2 . So 
from Eq. (|9|), the information per channel use in Eq. (|15p 
about 9 is 



MO) 



29 2 



(i - e 2 ) (> 2 + *b) 



(16) 



Our three channel identification schemes O, E, and B 
are just those considered in [H, H|. These schemes can be 
compared in different ways based on the quantum Fisher 
informations Jo(9), Je(9) and Jb(0)- The three schemes 
use different resources to different degrees, suggesting dif- 
ferent accountings. Schemes O and E involve just one 
channel use; scheme B requires two uses. Schemes E 
and B use entanglement, whereas scheme O does not. 
Schemes E and B involve a dimensional extension of the 
channel (with a more complex measurement); scheme O 
does not. Both Fujiwara [l[ and Sasaki et al. 0], com- 
pared these schemes on the basis of channel extension. 
Our Jq(9), J~e(9), and Jb(#) are information gained per 
channel use. Thus Fujiwara effectively compared 2Jq,(9), 
Je(#) and 2J-q(9). He found by his accounting that no 
one scheme uniformly dominated the others; each scheme 
was optimal for different 9. 

We compare our identification schemes O, E, and B 
for the general qudit depolarizing channel on the basis of 
quantum Fisher information per channel use. This com- 
parison emphasizes the essential cost of the channel to 
estimation. From a resource perspective this seems most 
natural, because if channel use is not relatively costly, we 
would always just probe the channel by scheme O with as 
many probes as needed to achieve the desired precision. 
Comparing quantum Fisher information per channel use, 
we find that 



MO) > Jo(0) and 
MO) > MO) 



(17a) 
(17b) 



for all 9 and all channel dimensions d, as shown in Fig. [5] 
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FIG. 2. Quantum Fisher information per channel use: Scheme 
E (upper solid line), Scheme O (dotted), Scheme B (dashed), 
Scheme 0,2 (solid lower). 



Comparing just schemes O and B, Fig. [3] shows that, 
above a d-dependent threshold for 9, scheme B always 
yields more information (per channel use) than does 
scheme O and that below 9 = l/\/3 scheme O always 
dominates scheme B. 

Ancilla depolarization: Identification scheme E uni- 
formly dominates schemes O and B for the qudit de- 
polarizing channel. Scheme E, though, may be unreal- 
istic in assuming that the ancilla qudit in the entangled 
pair undergoes no degradation. Suppose we model the 
evolution of the ancilla qudit by a depolarizing chan- 
nel T v with depolarizing probability 1 — ?/. If the probe 
and its ancilla qudit are initially in a maximally entan- 
gled state a = (see Eq. (fTTj) ). the channel output 
pe = {T v <g> T )(a) is 



Pe 



l-7]9 



r\9a. 



(18) 



This is identically the parametric family of states that 
would result from passing a Gp-dimensional qudit in a 
pure state a with no external entanglement through a 
depolarizing channel with parameterization h(9) = 7]9. 
Therefore, from the quantum Fisher information per 
use of Tg with this scheme E,, is 



Je v (0) = 



(19) 




Scheme O 



Scheme B 



i 



FIG. 3. Unentangled probes (scheme O) versus maximally 
entangled probe pairs (scheme B). Scheme B yields more in- 
formation above a dimension-dependent threshold. 



Unlike identification scheme E = Ei, scheme E^ for 
i] < 1 does not uniformly dominate schemes O and B. 
In fact, for any degree of depolarization r\ < 1 in the an- 
cilla system, there exist 9 where Je^(9) < Jo [9) and/or 
Je v (9) < Jb(0)- Simple calculation shows, for example, 
that Je„(0) < Jo{0) if 6 > g(rj) where 



9(v) 



ri 2 (d+l)-l 



?7(1 -v){d 2 -2) + if 



(20) 



Since g{rf) < 1 for any rj < 1 and any d, there is always 
s£ large enough such that Je„(^) < Jo{8)- Similar state- 
ments can be made regarding Je (9) relative to Jb(#). 
The broad point here is that, in this as in many other 
situations, realizing the theoretical advantage of entan- 
glement depends on precise control of the experimental 
apparatus, and particularly in this case, evolution of the 
ancilla qudit. 

Probe re-circulation : One might ask whether any ad- 
vantage can accrue to circulating the probe back through 
the channel, this being the scheme where the probe is 
passed through the channel, kept intact (no measure- 
ment), passed through the channel again, and, finally, 
measured. Conceivably, the probe could be re-circulated 
through the channel any number of times before measure- 
ment. Our schemes O, E, and B all accommodate probe 
re-circulation. The qudit depolarizing channel satisfies 
Tg n = Tgn where n is the number of channel passages so, 
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from Eq. ©, 

Jo,n(0) 

J E ,n(0) 



9 2(n-l) 



(W^) 

(1-0") ( " + ^t) 

2n ^2(2n-l) 
(1 - 02«) + _,L_ 



(21a) 
(21b) 
(21c) 



where Jo,n(0), Je,h(Q) and JB,n(#) are, respectively, the 
quantum Fisher informations per channel use (n channel 
uses with schemes O and E, 2n uses with B) about 9 
available by schemes O, E and B by measurement after n 
circulations of the probe through the channel. Some el- 
ementary calculus shows that Jo,n+i(0) ^ Jo,n(fi) f° r 
all n ^ 1, d ^ 2 and ^ 9 ^ 1, meaning that 
no degree of probe re-circulation is advantageous with 
scheme O. Since Je,ii(0; d) = Jo,n(#; d 2 ) and Jb,u(9; d) = 
Jo,2n(9',d 2 ), the same must be true for schemes E and 
B. The information J 0,2 (9) with scheme 0,2 is shown 
in Fig. 2 for various channel dimensions d. We see that 
with increasing dimension the entanglement schemes E 
and B yield information increasingly close to schemes 
0,1 and 0,2 without entanglement: Je(9) — > Jo i(9) 
and J B (0) ->■ J o . 2 (0). 

We observed that probe re-circulation, with any of our 
schemes, can only diminish the quantum Fisher infor- 
mation in the channel output. These observations were 
made, though, on the basis of quantum Fisher informa- 
tion per channel use. If channel use is not a significant 
cost, probe re-circulation can be advantageous. In fact, 
when the channel depolarizing probability 1 — 9 is small, 
re-circulating the probe just once can yield more infor- 
mation than entanglement. For the qubit depolarizing 
channel with 9 = .9, for example, the quantum Fisher 
information gained by measuring a once re-circulated 
scheme O probe is 2Jo,2(-9;2) = 9.42, while that gained 
by measuring a maximally entangled probe with no re- 
circulation is Je.i(-9;2) = 8.11. Entanglement clearly 
holds promise for channel probing, but depending on the 
relative cost of this resource, probe re-circulation or some 
other similarly simple classical strategy may be more ef- 
fective. In the case of the depolarizing channel, since 
probe re-circulation may be, for example, just a matter of 
maintaining the probe in the target depolarizing medium 
for some longer period of time, entanglement may need 
to be very inexpensive to be competitive. 



IV. PROBING WITH PARTIAL 
ENTANGLEMENT 

Various quantum Fisher informations J(9) in the out- 
put of the depolarizing channel probed by qudits with 
zero and maximal entanglement were presented and com- 
pared in the previous section. Here we consider qudit 



probes just partially entangled with an ancilla. The ef- 
fect of partial entanglement is not obvious; channels are 
known for which a partially entangled probe can yield 
more information than a fully entangled probe [|[. The 
essential conclusion of this section is that, for the qu- 
dit depolarizing channel of any dimension, partial entan- 
glement always yields information intermediate between 
Jo{9) from an unentanglcd qudit probe and Je(9) from 
a maximally entangled probe. 

The previous section's approach, represented by use of 
Eq. (|9|), is applicable to certain special kinds of partial 
entanglement — qudits which are fully entangled on just a 
subspace. We briefly pursue this point and then turn to 
the question of partial entanglement generally. Consider 
a pair of qudits in the partially entangled pure state a = 
I \T/) ( \I/ 1 with Schmidt decomposition 



I*) 



V>i|n) + ... + ^|rfd) 



(22) 



where tpj ^ and tpf H h = 1. Suppose that d of 

the ipi are \j \fd~ and that the remainder are zero, where 
2 < d Q < d. These qudits are maximally entangled on 
a subspace of dimension d < d with no entanglement 
outside this subspace. The output of the depolarizing 
channel probed by the first of these qudits is 



1 



dd D 



■I d ® l do + 9a 



(23) 



and associated with this output is the quantum Fisher 
information 



J{8) = 



1 



(24) 



This shows for qudits with this special form of partial 
entanglement that 



Jo(9) < J{9) Me). 



(25) 



We now prove that Eq. (|25|) is true generally for any 
partial entanglement. 

Suppose two qudits are partially entangled generally 
as in Eq. ([22]) . Probing the channel with the first qudit 
produces the channel output 



Pe = (r e ®I)((7) 



1 - 9 



I d ® D + 9a (26) 



where 



(27) 



operates on a single qudit. Redefining D := 
Si=i V'iNK**! as an operator on both qudits gives pg = 
Pc + Pq where 

1 



PQ 



D + 0|*)(*| and p c 



-C (28) 
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are positive scmi-dcfinitc operators with orthogonal sup- 
ports S — {\ii),i = 1, ...,d} and S = {\ij),i 7^ j}, 
respectively, determined by the Schmidt decomposition 

Eq. (HlD,andC = E^^ 2 Kj>(ul- 

A general result, whose proof is straightforward, is 
that, if pg = pw + p2g where p\g and p2S have orthogo- 
nal supports, then the score operator associated with pg, 
i.e. L, can be expressed as L = Li + L2 where Lj is the 
score operator associated with pig. Also the support of 
hi is orthogonal to that of Lj and pjg if i ^ j. Thus 
L = Lq + Lc where \jqPc = LcPq = such that Lq 
and Lc satisfy, respectively, 



L QPQ + p Q L Q = 2dg PQ = 2 ( |*)(*|--D ) , (29) 



Lcpc + Pc^c = Zdgpc = - jC. 

a 



We have immediately from Eq. (|30|) . then, that 

1 _ 



1 



(30) 



(31) 



where Ic is the identity on S . While no similarly simple 
solution exists for Eq. (|29[) . substituting pq from Eq. 
in Eq. (J29J) yields 



1 1 



(10 



dO 



Tr(DL Q ). 



(32) 



We now derive, using our results from Eqs. (|29l 
and (|30p . an expression for the quantum Fisher infor 
mation J(0) in the channel output Eq. (f26|). 

J(0) = Tr(Ld e p e ) 

= Tr (Lcdgpc) + Tr (L Q dgp Q ) 



Tr 


L c ( 


d 


- 1 


d(l 


-f) 


d 


- 1 


d(l 


-9) 



Tr 



(*|LqI*> 

d- 1 1 



1 



*)(*|--D 



dO 

Tt(DLq) 

- 1 1 
~ dO 



dO 



Tr(DLq) 
Tr (DL Q ) 



d0{\ - 0) 
d-l 



Tr(DL Q ). 
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l T Dx 

d0(l - 0) dO 

where x is the d dimensional vector 



( L A 

L22 



(33) 



(34) 



\Ldd) 



with La = (m|Lq|m), D is, now, a d x d diagonal matrix 
with diagonal entries ipf , . . . , if}^ and 1 is the d dimen- 
sional column vector of ones. Expression (|33[) shows, sig- 
nificantly, that only Lq's diagonal elements are needed 



for J(0). In appendix [A] we demonstrate that 

T d-l d 2 

l T Dx = 



1 - 



dO 



(1 - o)(i - e + doy 



b T z- x b 



(35) 

where b is a d(d — l)/2 dimensional vector with lexico- 
graphically ordered entries 

/ ipiip 2 \ 
b = 2 (36) 

\lpd-XtpclJ 

and Z is the d(d — l)/2 x d(d — l)/2 matrix 

Z = J + Q$7 T r2 (37) 

where a = d6/(l — + dO). The two constituent matrices 
of Z are 



Vf + V'l 



\ 



V 







(38) 



and a d x d(d — l)/2 matrix = fd{i^i, 
fd(') is defined recursively by 



fd+\{x,T d ) = 



and 



r d 

xl d fd(*d) 



. . , ijJd) with 

(39a) 
(39b) 



where r<j = (n . . . r d ) and is a row vector of (d — l)(d- 
2)/2 zeros. For example 



n 4 = 74(^1,^2,^3,^4) 

M 2 ^ 3 V 4 

V>i ^3 ^4 

Vi ^4 

\0 ^ 2 Vs> 



(40) 



Combining Eqs. (|33| and (|35|) yields 
1 



J(0) 



(i-*)(* + a^i) 

d 

+ (i - 6>)(i -o + doy 



(41) 



While Z is singular in cases where the qudits' entangle- 
ment is confined to a proper subspace, b T Z _1 b is always 
well-defined. With this in mind one can readily check 
that, for maximal entanglement on a subspace, Eq. (|4ip 
is identically Eq. flUj) , 

To prove Eq. (f25j) we must prove that 



b T Z _1 b < (d- 1)- 



1 - 



d0 



d 2 e» 



(42) 
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or, equivalently, that 



b T (an T fi + J) _i b < 



(43) 



d-1 



Let A m in[M] and A max [M] denote, respectively, the min- 
imum and the maximum eigenvalues of a positive semi- 
definite matrix M, and let denote the d x d matrix 
Arf = r2J _1 f2 T . Let ip be the length-d column vector 
with components ipj. U sing the relation b = n T ^, the 
matrix inversion lemma [22[, and the Rayleigh-Ritz the- 
orem [231, we have 



b T (afl J fl 



j) : b = i> T (i( a n T n + j)- 1 ^^ 



{al d 
al d + 



max 



a + A min [A rf ] 



a + A m ax 



[Ad] 



(44) 



According to Eqs. (f4"3"|) and (JUJ), we only need to show 
that A max [Ad] = d— 1 to prove the inequality of Eq. (|2"5f . 
First we show that d — 1 is an eigenvalue of A^; then we 
show that it is the maximum eigenvalue. Let A be an 
eigenvector of A^ with associated eigenvalue A so that 



A d A = AA. 



Let Ai = J~ L fl 1 A. Then 



j- 1 n T nA 1 = AA! 



(45) 



(46) 



so A is also an eigenvalue of J _1 fi T ri. One can check by 
direct substitution that Eq. (|46|) is satisfied by A = d — 1 
with (unnormalized) eigenvector 



Ai 



f 



1 



■01^2 1pl1p3 i>d-\^d 



(47) 



In appendix [B] we prove that d — 1 is the maximum eigen- 
value of Ad, and this completes the proof of Eq. (|25l) . 



V. OPTIMAL MEASUREMENTS 

Realizing the precision promised by the quantum 
Fisher information in the Cramer- Rao inequality depends 
on optimal measurement of the channel output. The 
quantum measurement that saturates the inequality is 
not always obvious, and in some cases — certain vectors 
of parameters, for example — no such measurement exists 
jig ]. For the qudit depolarizing channel, is is possible to 
describe a measurement which yields a (classical) Fisher 
information which equals the quantum Fisher informa- 
tion. For this channel the density operator after evolu- 
tion is such that its eigenstates are independent of the 



channel parameter (only the eigenvalues depend on the 
channel parameter), i.e. 



= 5>(0)hX*il 



(48) 



where Xi(9) and \sj) are the eigenvalues and eigenstates 
respectively. Then [18| the quantum Fisher information 
is 



(49) 



However [18J , this is exactly the expression for the classi- 
cal Fisher information for a probability distribution with 
probabilities {K{0)} . Thus the quantum Fisher informa- 
tion can be attained via a measurement which results 
in this probability distribution. Clearly measurement in 
the eigenbasis of pg accomplishes this. Here, one of the 
eigenstates of pg is the initial state \p) of the pair of qu- 
dits prior to action by the depolarizing channel. It follows 
that the projective measurement for which the projectors 
are 



n = |/i)(/x| and 

n x = j- |/*></*| 



(50) 
(51) 



yields a probability distribution, which gives a classical 
Fisher information exactly equal to that of the quantum 
Fisher information. For an entangled input state such as 
that used in scheme E, this optimal measurement requires 
measurements involving entangled states. It is unclear 
whether there is a scheme which involves unentangled 
states at the measurement phase and which attains the 
quantum Fisher information; such a scheme exists for 
unitary parameter estimation [23j . 



VI. SUMMARY 

Wc have given an extensive treatment of channel iden- 
tification for the general qudit depolarizing channel, com- 
paring different probe preparation schemes based on the 
quantum Fisher information attainable by channel prob- 
ing. We showed for depolarizing channels of any dimen- 
sion that maximally entangling the channel probe with 
an external qubit uniformly realizes the maximum theo- 
retical advantage, but that this advantage disappears if 
the ancilla qubit is not sufficiently shielded from depolar- 
ization. Wc primarily compared the same three channel 
identification schemes considered in Q, though our 
results apply generally for dimension d < oo. We de- 
parted from |l|, l2[ most significantly in choosing to com- 
pare identification schemes on an information per chan- 
nel use basis. This is a practically meaningful basis of 
comparison; in any case, our results for each scheme are 
readily amenable to other comparisons. 

Wc have considered here identification schemes involv- 
ing bipartite probe entanglement, with an ancilla qudit 



8 



or with another probe. Multipartite probe entanglement 
of various orders and types remains for study. We hope to 
be able to address at least some important multipartite 
entanglement schemes with our present methods. 



Appendix A: Diagonal elements of the score 
operator 



Substituting from Eq. 



1. 



l*X*l 



1 



2d 



into Eq. (gSJ) gives 
(LqD + DLq) 

(Lq|*)(*| + |*)(*|Lq).(A1) 



Acting with this on \jj) and acting on the result with 
gives 



1 



fe=i 

where = (m|Lq|jj). Note that the score operator 
is Hermitian and in this case real; thus Ljj = Lji. Thus 
Eq. (|A2|) is a system of d(d+l)/2 distinct linear equations 
for the components of Lq. Our aim is to invert these 
for the diagonal elements of Lq. Using x as defined in 
Eq. (J3U) and 

/ il2 \ 

L13 



Lid 
\Ld-idJ 



(A3) 



gives 



R S 
S T T 



(A4) 



where a is a d dimensional vector given by 



d-l 



(A5) 



b is given in Eq. (f3l))) , and the concatenation (x yj 
is the column vector consisting of the components of x 
followed by those of y. The submatrices R and S within 
Eq. (|A4[) are readily derived from Eq. (|A2[I and are 



R 



(pn 

P22 

\o 







Pdd/ 



1 - 



d6 



D 



(A6) 



and 



S = (9D 1/2 O n 



(A7) 



where fid is as defined via Eq. (|39b[) . The final submatrix 
is 



1 - 



(A8) 



and this can be verified by considering Eq. (|A2p for i < j. 
The left side yields the row of b corresponding to the 
indices % and j and the right side yields the following 
contribution from Ty 



1 



d d 

ipjipkLik + ^2 ip t il> k L jk I . (A9) 



The first term clearly corresponds to (1 — 9)/d Jy. To 
check the second term, we can show by an iterative pro- 
cess that 



fidy = 



(A10) 



Jlk^d ^k^dky 

Using this, we can again show iteratively that 



( E^i ^kLik + Efe# 2 i>kL2k \ 
E^i i>kLik + E fe ^3 ^kL 3k 

VEfc^-i i>kL d -ik + E^d ipkLdkJ 



(All) 



and this corresponds to the second term. This justifies 
the form of Eq. ifAi]) . 

Applying the matrix inversion lemma (22j to the solu- 
tion of Eq. (|A4j) for x gives 



(R - ST S 



1 C T\ 



(a- ST" 1 !}) 



= [R- 1 + R- X S(T - S T R- 1 S)- 1 S T R" 1 ] 
(a-ST _1 b) 

From (|AT2|) . and using 



R a 



d-l 



i-e + de' 



(A12) 



(A13) 



a lengthy calculation yields 
d 



Dx = 



i-e + de 
d 

i-e + de 
d-i 



(a - ST _1 b) 
S(T- S T R- 1 S)- 1 S T 



1 - 6 + dO ' 



1 - R^ST^b 



(A14) 
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Using the relations l T a = (d — l)/d and S T 1 = 6h, we 
calculate from Eq. (|A14jl that 



l'Dx = 



d- 1 



dO 



Recalling the definitions of S and 
and (IA8I). we calculate further that 



i-e + de {l-e + dey 

xb T (S T R _1 S — T) _1 b. (A15) 
T in Eqs. ([AT]) 



l T Dx = 



d- 1 



1 - + <20 



d 2 



(l-0)(l-6» + rf6») ; 



where 



Z = J 



ir2 T r2 



b'Z^b (A16) 



(A17) 



with a = d6/(l-6 + d6). This gives Eq. ([35]) 



Appendix B: Maximum eigenvalue of Ad 

The matrix = [a,-,-] has elements 



Ed 
k=l 



a + n 



if i 



3, 



Hi ^ j. 



(Bl) 



For example, 



V3 



•02+</>3 



</'5 



(B2) 



Define the angles 0^ for 1 $C i < j : ^ d by cos </>,j = 
V'i/yV'f + ^1 an ^ let be the length-rf(d - l)/2 vec- 
tor $ d = (</> 12 0i3 ... (j>d-\,d)- The matrix A d 
is a function A ( j(<i> ( i) of the angles in Consider 
the characteristic function h(X) = det(AI f / — Ad) of 
Ad- The eigenvalue A max [Ad] = d — 1 if and only 
if h'(X) is positive for A > d — 1. This is a conse- 
quence of the following. If the eigenvalues of A^ are 
A m i„, . . . , A max then h(X) = (A - A min ) ... (A - A max ) and 
h'(X) = /i(A)[l/(A-A mm ) + ...l/(A-A max )]. Clearly 
if A > A max then h(X) > and h'(X) > 0. A graphical 
representation of the characteristic function reveals that 
at some point between A max and the eigenvalue imme- 
diately smaller than this, h'(X) must be negative since 
ft. (A) in this region is between two of its roots. Thus if 
h' (A) > for all A > A then A = A max . We have 

d 

ft'(A)=^det(AI (i _ 1 -A^) (B3) 

i=l 

where A" is the matrix Ad with its ith row and column 
deleted. The summand in Eq. (|B3fl is the characteristic 



function of A", so we need just show that no A^ has 
an eigenvalue greater than d— 1. A standard calculation 
yields 

A^Ai-i^'j + A 4 (B4) 

where A* is a diagonal matrix with diagonal elements 
cos 2 c/>ij , and & d is the vector <J><j with the d — 1 angles 
4>ij and <j)ji deleted. Therefore, by Weyl's inequality [22| . 

AmaxtAjf] < A max [4i-l(^)] + A max [A*] 

+ (B5) 



We have A max [A 2 ($ 2 )] = A max [A 2 ($ l 3 )] = 1 by direct 
calculation from Eq. (|BT]) . Thus A max [A 3 l ] ^ 2 by 
Eq. (|B5jl . and A max [A3] ^ 2 by consideration of Eq. (|B3jl . 
Arguing recursively, we conclude that A max [Ad] ^ d — 1 
for d ^ 2. We saw that d — 1 is an eigenvalue of A^ so 
A roax [A d ] =d-l. 
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